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SYNOPSIS 


la this V/Ork angular distributions a'ld excitation cuirves 
of (d,j') reaction on C and Be have been obtained at the 
sab-GOUloinb bonbardmg energy in the range of 500 to 1800 KeV 

using the 2 MeV 7an-de-Graaf f accelerator at I I.T., Fanpur. 

12 

Strong resonances were observed in the case of C and compa“ 
ratively weaker resonances in the case of ^Be in the excita 
tion curves. The angular distributions in both cases /however, 
showed a peak at about 35^j^ typical of direct stripping mode. 

At resonant energies the main effect was in the enhancement 
of the backward peak. Some data on both these reactions was 
previously available. However, the analysis was done either 
in terms of compound nucleus process alone or the stripping 
mode alone. In the case of compound nucleus mechanism, it was 
necessary to assume contribution from rather high partial waves, 
and in the case of stripping mechanism, it \7as necessary to 
assume optical model parameters of DVJBA calculations which 
were not consistent with those used at higher energies where 
the angular distributions shov/ better strapping patterns. It 
was, therefore, decided to analyse the data in terms of coherent 
contribution of the reasonant and non-resonant amplitudes and 
then compare the non-resonant cross-section wa th the DWBA cal- 
culations with optical model paraineters consistent with those 



used at higher energies For a proper consideration of the 
interference effects it was necessary to laj e the data at 
recfuired energies and angles/ although soi c oata was availa- 
ble from previous work 

The amplitude was assumed to be of the lorm 
^ _ ^Resonant ^ ^Non-resonant 

in the vicinity of resonance pea] Where j^-f^®sonant ^^j^^^^luded con- 
tribution from the main peak and one or two neighbouring peaks 
It was assumed that the energy dependence of the resonant 
amplitudes was given by Breit-Wigner terri/ with phases of the 
neighbours relative to the main resonance The direct ampli- 
tude was considered relatively insensitive to energy over the 
region of the resonance and was assumed to be at the most 
linearly dependent on it with the undetermined phase relative 
to the main resonance This assumption about the total ampli- 
tude leads to an exnressioh for the total cross-section with a 
term corresponding to the non-resonant cross-section, and 
others depending upon energy in the specified way The expe- 
rimental cross-section in the neighbourhood of resonances in 
the excitation curve Was then subjected to a least squares fit 
using the expression for cross-section with strength^ of various 
terms reso’^ant energies and widths > as parameters The feast 
squares fit prov-des a -^lue of the non-resonant cross-section 
at each energy of- the excitation curve obtained at each angle 



It IS then possible to reconstruct tne angular distribution 
due to the non-resonant process alone 

These extracted non-resonant angular distributions were 
compared with a D! BA calculations of the stripping process 
made with a program DWUCK available at the Dhabha Atomic 
Research Centre It was foxind that there were satisfactory 
agreement with the shape and absolute cross-section between 
the two 

Ih IS/ therefore concluded that while the contributions 

of various processes in determining the cross-section are to 

be determined by coupled chanral calculations the simpler 

model of adding resonant and non-resona^t amplitudes suggested 

physical considerations in the classical limit seems to 

yield satisfactory results Uith this approach the contribu- 

12 

tion of the direct reaction in the case of C was found to be 

o 

about 20?^/ and that in the case of ■'Be about 90/ In view of 
this, the assignment of spins and panties made previously by 
neglecting the non-resonant part needs to be re-examined An 
attempt to fit the extracted direct cross-section with Coulomb 
distortions alone did not provide satisfactory fits inspite of 
the low bombarding energy Nuclear distortion is therefore, 
found important even at these energies Ihe finite range and 
non-local corrections in the DWBA calculations did not produce 
any major ch^ge in the shape However, the absolute cross>- 
sectJon reduced by about 10% with these corredtions applied 
to 'the zero range apprtecimation 



CHAPTER I 

INTRODUCTION 

The striking feature that was noted in the early stages 
of nuclear reaction work was the appearance of narrow resonances 
in the cross-section as a function of bombarding energy* This 
led to the suggestion of the compound -nucleus mode of nuclear 
reactions by N Bohr/ ^ with compound-nucleus life time of the 
order of 10 lo”^^ sec which was much larger than the typi- 

cal transit time of the bombarding particle across the nuclear 
diameter The theory of nuclear reactions on the basis of com- 
pound-nucleus formation mechanism was siabs equently developed 
by many authors and was successfully applied to determine the 
properties of many compound nucleus levels 

In the case of nuclear reactions proceeding via compound 
nucleus formation mechanism xt is expected that excitation fun- 
ction will show resonances "Hbose shapes are given by the 
Wigner forinUla Since the formation add decay of the compoxnd- 
hucleus are independent of each Other the same reaction jprod-uots 
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would be obtained even if a particular level of the compound- 
nucleus IS formed by two different initial channels The angu- 
lar distribution of the reaction products proceeding through a 
single level is expected to be symmetric about 90° in the centre 
of mass syston Since the shape of the angular distribution is 
characteristic of the partial waves in the outgoing channel, it 
IS possible to make an inference about the spin of the compound- 
nucletis level If at relevant excitation energy two or more 
levels of the compound nucleus of different panties contribute, 
the angular distributions are not expected to be symmetric 
and the shape of the angular distribution is expected to vary 
noticeably with the bombarding energy The width and the den- 
sity of the compound-nucleus level are expected to increase with 
excitation energy and when a large number of levels contribute 
it IS expected that, due to essentially random phase of contribu- 
tion the angular distribution would again be symmetric and 
weakly dependent on angle 

Even before the suggestion of compound -nucleus process 
by Bohr^, Oppenheimer and Phillips^ had proposed a radically 
different model to account for the higher probability of (d,p) 
reactions as compared to (d, n) reactions They pointed out 
that at the low bombarding energy which was used at that time,> 
tb® proton in the deuteron wduM not cross the Coulomb-barrier of 
tne target whj.le the neutron would get close enough to get cap«- 
tub^- It was thus possible to ^fplaln the preferential proton 
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emission which is just the opposite of the predictions of com- 
pound-nucleus theory The full implication of this model was 
not clear at that time since at very low deuteron energies the 
angular distributions were isotropic as only s-wave interactions 
took place In the fifties/ (d p) and (d/ n) reactions were 
studied at higher energies and they were found to be equally 
probable, showing that coulomb force pjays the secondary role 
However, the most interesting feature was the angular distribu- 
tion The cross-section peaked either in the forward direction 
or at some small angle and decreased to small values at back- 
ward angles These features were independent of bombarding ener- 
gies, and peak position depended upon the target used These 

4 

features were explained first by Berber and later more accura- 

5 

tely by Butler The explanations were based on the single step 
model of the reaction in which the deuteron is stripped of one 
of its nucleons, due to the low binding energy of the deuteron, 
during the close passage of the captured nucleon relative to the 
tafget nucleus The other nucleon can still have a relatively 
large impact parameter in view of the large deuteron radius 
in the stripping process therefore the deuteron as a whole and 
the target do not combine to form a compound nucleus Butler 
showed that the angular distribution is essentially of the form 
o£ the square of the spherical Bessel fianction ^ Of order 

1 ^^^-where 1 is the orbital angular momentum transfered to the 

nucleus’' by the captured nucleon Q is the linOar momfentura 
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transferred and R is the nuclear radius The study of deuteron 
stripping reactions therefore provided a new technique of nuclear 
spectroscopy to study the properties of the low lying levels of 
the residual nucleus, while the reactions proceeding through 
compound nucleus mechanisms led to the study of highly excited 
levels of the compound nucleus The inverse (p, d) and (n, d) 
reactions are called the pick-up reactions in which the inci- 
dent nucleon picks up an appropriate peripheral nucleon to form 
the emergent deuteron without fo3nning the compound nucleus The 
Bui ter theory is directly applicable to such reactions and the 
peak position in the angular distribution now depends upon the 
angular momentiim of the picked nucleon in the target nucleus 
The pick-up reactions thus provide information about the ground 
state of the target nucleus 

In a rearrangement collision the amplitude of the reaction 
7 

a P is given by 




= / e Vg f dT 


(1 1 ) 


where is the internal wave function of the reaction products 


in channel g, the final state interaction and, Vthe total 
wave function In the Cd,p) reaction Vg = where 

is the pr^ ^y -neutron interaction and is th#' sum of 
the interacti^^ of the proton with the nucleons of the target 
•^e pure stripping model it is assumed that the interaction 
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8a 

can be negletced. The Butler treatment and the other 

Sb 

alternative treatments that appeared soon afterwards are essen- 
tially equivalent to making a plane wave Bom approximation 
(PWBA) in which the total wave function '1*13 replaced by plane 

wave e *’ a m the incident channel. The method used by 
8a. 

Butler involves v;ave function matching at the nuclear boundary 
and it leads” to same result as the PWBA with the difference that 
the radial integral for amplitude is cut-off at the nuclear radi- 
us R and does not extend to the origin. 

In spite of the success of the P5*7BA theory for deuteron 
stripping reactions, the theory does not yield a particularly 
detailed account of the angular distributions For medium and 
heavy nuclei the agreement between the PWBA theory and the expe- 

g 

rimental angular distributions is poor. Even for light nuclei 
there are cases in which the assignment of the orbital angular 
momentiam is ambiguous. Further, PWBA gives a value of the 
absolute cross-section which is several times the measured cross- 
section. 

The exact expression (1.1) for the amplitude of the 
reaction can also be put in the form 

= /Xg”'\rp)fp (Vp -Vg) y dr ...(1.2) 

where is a model potential and Xg(fg ) , represents the relative 
motion in the model potential. The PWBA is improved upon 
by using the Distorted Wave Bom app:^imation (DWBA) theory 
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The total wave function ¥ in expression (1.2) is replaced 
not by the plane unscattered wave but by the scattered wave 
X„ the incident channel. The scattering state wave funC' 

tion is generated from the optical model, where the optical 
model parameters are determined from the available scattering 
data in the incident channel. The model potential Vg is chosen 
equal to the optical potential which gives correct scattering 
in channel B . Thus ^ scattering-state wave function 

in channel 6 . This wave function is generated from the optical 
model potentials which determine the optical model parameters 
for a given system. As mentioned before, for the (d, p) reac- 
tion VftSsV^ + V therefore, 
p pn pN 

= V ‘V" ...(1.3) 


It IS assumed that Vgis a close approximation of , 

which IS the total interaction the proton experiences due to 
the nucleons in the target nucleus. Thus in the DWBA theory the 
plane wave is replaced by the distorted wave, which is a better 
approximation to the total wave function and instead of neglec- 
ting the interaction of the outgoing nucleon with the target, it 
IS assumed that is equal to the optical interaction to a good 

approximation . 


The optical model potential, which is in general required 
in DWBA treatment, includes both Coulomb and nuclear interactions 
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The inclusion of nuclear interaction in both the incoming and 

outgoing channel represents a departure from the pure stripping 

mechanism as envisaged originally by Butler. It implies that tine 

deuteron as a whole does experience nuclear potential. Similarly 

the outgoing proton does not entirely pass through the region 

exterior to the nucleus. However, the optical potential is 

adjusted to give the correct energy -averaged scattering 

amplitude and does not reproduce the fluctuations in the elas- 
12 — 

tic channels. The potential Vg is not identical with 
DWBA treatment, therefore, does not take into account the possi- 
ble compound nucleus formation and subsequent contribution to 
the reaction leading to narrow resonances in the cross-section. 
DWBA, however, provides a more accurate description of the dyna- 
mical behaviour of the incoming deuteron and outgoing proton 
and by doing this provides a sounder basis for using proton- 
neutron interaction as the operator for the matrix elanent in the 
amplitude (1.2). 

Generally speaking DWBA calculations require extensive 

use of computers. Several computer programmes for such calcu- 
. 13 

lations are now available. Many programmes use a zero -range 
approximatign for the neutron-proton interaction since it re- 
duces the dimensionality of integral (1 .2 ) . : CalGulations 
which take into account the finite-range of are lengthier. 

Programmes which include such finite range dorr ecti on as well 
as the spin-orbit term in the optical potential are also 
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available. Computations are also made witb non-local poten“ 
tials. There is usually a provision to use a cut-off radius 
in integral (1.2) as an option. The DWBA calculations usually 
give a much more satisfactory account in the entire range of 
angular distribution as compared to PWBA calculations. Further, 
the absolute cross-sections predicted on the basis of DWBA cal- 
culations are in far better agreenent with the experimental 
results . 

For bombarding energies which are lower than the Coulomb 
barrier the (d,p) cross-sections are expected to be small and 
the Coulomb distortion is expected to be the predominant effect^"^ 
Secondly/ the relative contribution to the reaction due to compound 
nucleus formation is expected to be larger This is due to the 
fact that the relative strength of the low partial waves in the 
incident wave increases as energy is decreased and the compound 
nucleus formation is more probable with such low partial waves 
with correspondingly low centrifugal barrier On the other hand, 
when compound nucleus is formed at high bombarding energy it is 
highly excited and therefore has several open channels. The 
fractional contribution of the decaying compound nucleus in the 
outgoing proton channel is therefore reduced. At low bombarding 
energies / therefore, it is expected that the cross-section would 
be low, distortion effects will be primarily due to Coulomb field 
and the compound-nucleus process will malce a significant contri- 
bution to (d,p) reaction. 
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The Q value of (d,p) reaction is an important conside- 
ration particularly at low-oombarding energies When the Q 
value of the reaction is large the momentum of the outgoing 
proton is large. To have this large momentum it is expected that 
the proton is close to the neutron in the deuteron when the 
neutron is stripped. Therefore, such protons pass close to the 
nuclear region and experience marked nuclear distortion. It is 
expected, therefore, that in the high Q reaction the plane wave 
approximation will be poor. This has been observed for several 
cases. Tne (d,p) reaction on and ®Be have been 

studied at bombarding energies in range 0.8- 2.0 MeV. In these 
reactions Q values are high. Therefore, both compound -nucleus 
effects and high distortion effects in (d,p) reactions are 
expected. 

When both compound -nucleus and direct amplitudes are 

present, it is expected that they will add coherently. While 

-21 

the duration of a direct event should be of the order of lo 

4 

seconds and that of the compound -nucleus some 10 times greater; 
the duration of the wave packet which is incident from the acce- 
lerator IS of the order of lO seconds. The difference in 
delay between the direct and compound-nucleus mechanism is 
therefore insignificant. The well known experimental evidence 
of the existence of interference is that between the potential 
scattering and resonant scattering of neutrons . The time 
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intervals involved in direct interactions are the same as those 
involved in potential scattering. VJhen the interfering ampli- 
tudes are fe\7 in number the interference effects will be noti- 
ceable in the cross-section measurements. VJhen the number of 
interfering amplitudes is large the interference terms may tend 
to cancel. 

In the presence of an isolated resonance in the (d,p) 

on n 

reaction one may assume the presence of two amplitudes q 
and q adding coherently. If the resonance is not 

well isolated one may have to add the amplitudes due to neigh- 
bouring resonances. In the statistical region where a large 
number of resonances contribute, this approach is not possible. 

A 

In the analysis of excitation function in the presence of one 

or more resonances, the resonant amplitudes may be assumed to 

2e 

vary according to the Breit-Wigner formula while the Non- 
Resonant (nr) amplitude may be assumed essentially constant 
over the width of the resonance. It is then possible to extract 
the NR cross-section from the experimental excitation curves by 
fitting them with the algebraic function written down on the 
basis of above assumptions. The extracted NR cross-section can 
then be analysed on the basis of DVJBA. 

In the present work we have studied deuteron induced 
reaction at low bombarding energies on light targets with, the 
2MeVVan de Graaff accelerator facility available at the Indian 
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Institute of Technology, Kanpur. and ^Be were used as 

targets Previous data on carbon-12 (d,p) reaction in this 

17a 

energy range have been published in 1960 which showed pre- 
sence of resonances and the data was anlysed by the authors 
in the frame-work of compound -nucleus process alone, although 
the angular distributions did retain some feature of the mam 
stripping peak which shows up more prominently at 4 
A preliminary calculation of pure stripping at this bombarding 
energy showed that the stripping cross-section may not be negli- 
ble. It was, therefore, decided to study the contribution of 
the direct mode relative to the compound-nucleus mode, by consi- 
dering coherent addition of amplitudes. For this, it was nece- 
ssary to check previous data and to take additional data to 

study the interference effects. Some data on (d,p) reaction on 

9 1 9 

Be in this energy range was also available which was analysed 

by the authors in terms of compound -nucleus process alone or 

20 

the direct reaction process alone using optical parameters 
not consistent with those used at higher energies. 

12 

In the present v'ork it was found that in C(d,p) reac- 
tion the direct mechanism makes about 20% contribution as compa- 
red to the figure of less than 2% mentioned in the earlier 
1 Tlo 

worK The extracted NR cross-sections show strong distor- 
tions as expected in relatively high Q value reactions. 

The extracted cross-section can not be fitted by pure Coulomb 
distortions. Reasonably good fits are obtained with optical 
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potentials with parameters consistent with those used at 
higher energies In the case of beryllium-9 the compound- 
nucleus contribution is found to be small. The extracted 
cross-sections can again be fitted reasonably well with opti- 
cal model parameters consistent with those used at higher 
energies « 



chapter II 

NUCLEAR REACTIONS 


Introduction 

In a nuclear reaction a target nucleus 'A’ with n^ 
number of nuclei is bombarded with a projectile 'a' consis- 
ting of n^^ number of nuclei and the reaction products are 
detected in the asymptotic region. The reaction 

a + A 5^" B+b ...(21) 

IS described by denoting the initial and final channels by 
a and 6 respectively. Besides the identification of initial 
and final products other relevant quantum numbers are denoted 
collectively by symbols i and j respectively, so that the reac- 
tion IS represented by (a,i) — ► (^/j) when necessary. 

The total system thus consists of n,+ n^= N particles 

with the Hamiltonian H and a vzave function 'i' - In practice the 

state IS a non- stationary state, to be represented by colli- 

21 

Sion of wave packets, but it is generally represented as a 
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stationary state satisfying 

H ¥ = E 'F ... (2.2) 

The stationary state wave- function '? can he expanded 
in terms of a suitable basis such as independent particle 
model (IPM) states of N particles- It may be noted that for 
sufficiently high excitation energies, the finite well JPM pro- 
vides both bound and unbound stafes- The excitation energy 
can either be distributed over many excited single particle 
bound states or it may be in the form of one or few particles 
excited into the single particle continuum. The bound states 
of IPM therefore provides a discrete spectrum in the same 
region where the unbound states provide continuum. The unbound 
states represent the flux of particles in the asymptotic region. 
In nuclear reactions the energy of the total systan is suffi- 
ciently high to provide unbound components at least in the 
initial channel. The wave- function ¥ could then be written in 
the form 


bound unbound 


.. (2.3) 


where ’F, is a summation over the bound basis states and 

bound 

m ^ IS a summation over unbound basis states, 

'F unbound unbound 

is often expressed as 


'F 

unbound 



..(2.4) 
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where is the product <t>} 3 (j^of the internal wave functions of 

clusters b and P and ^ (^o) describes their relative notion. 

B P 

While this description of could be used for large 

values of r^ , it is sometimes used for small values of the 
order of nuclear radius as a model wave function- It is assumed 
that in the asymptotic region where the bound states do not 
contribute the total wave function has the form 


"" —I 


B 




3 2TTh 


2 % ^a3 


. (2-5) 


This expression describes an incident plane wave in 
channel ct and outgoing spnerical waves in various channels 3 . 
Since only the unbound states contribute in the asymptotic 
region/ the expansion (2.4) must have the asymptotic form (2»5). 
The quantity (i^/Fg) is called the amplitude of the reaction 
a Band the cross-section is proportional to the square of 
this amplitude. Measurement of cross-section therefore provi- 
des a sampling of the unbound component of the total wave- 
function . 

’^bound '^unbound' however, are two components 

of the total wave function and these are directly coupled with 
each other through the Schrbdinger Eq. (2-2). Thus, if we 
introduce projection operators P and Q with the property 

= Vbound ' Wd ...(2.6) 
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the Schrodinger equation directly leads to the coupled 
equations : 

( P \ =r PT V 

PF'^ unbound PQ bound 


'^'bound 


H V 

QP unbound 


. * . (2.7) 


where 

Hpp s PHP s QHQ 

= PHQ Hqp s qhp *..(2.8) 


'^bound ^ superposition of bound IPM states and the 
strength of such components are expected to vary sensitively 
with the system energy leading in turn to energy sensitivity 

''bound- cue to the coupling of Vbound ''bound- 

energy sensitivity will also appear in therefore 

in the reaction amplitude. If the bound component of the 
total wave- function could be neglected, then the '^'unbound 
expected to be relatively energy insensitive since they re- 
present continuum states of the system. 

In the nuclear reaction theory exact calculations are 
at present not attempted, and it becomes necessary to use appro- 
ximations and models. In the direct reaction theories the bound 
component of v/ave function is not c'onsidered explicitly. The 
form (2. 4) for the unbound wave function is usually used as 
the total model wave ; function in the entire region including \ 
the nuclear interior. The model potential Vg which determines ^ 
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?g(rg) IS empirically determined to fit the observed data in 
the corresponding scattering channel alone. It is in general 
not possible to find model optical potentials which will re- 
produce narrow fluctuations typical of a system with many 
degrees of freedom. 


The approach of the compound nucleus theory is to assume 
that^ while ¥ , ^ is the dominant term ofF in the region 

starting just outside the nuclear volume, i® the domi- 

nant term inside the nuclear volume. The coupling between the 
two IS then taken into account by the physical requirement of 
matching the two wave functions at the nuclear surface. The 
fluctuations can then be understood in terms of rapidly chan- 
ging boundary conditions as determined by assump- 

tion of smooth behaviour of the logarithmic derivative on the 
nuclear surface then leads to Breit-Wigner formula for the 
resonance The resonance is attributed to a predominant exci- 
tation of a bound state component. 


Both the treatments involve approximations to the actual 
problem, end introduction of models within these approxima- 
tions. The support for introduction of such simplification 
however^ is provided by the experimental data which in diffe- 
rent suitable circumstances correspond to good approximations 
to the two different treatments. Further, the physical consi- 
derations at least in the classical limit lend support to the 
introduction of these simplifications. 
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In a semiclassical approach the nuclear reaction can 
be considered in terms of the development of rearrangement 
collision in time. Thus, when a deuteron is incident on the 
target there is a possibility of the neutron being captured 
and the proton released. Alternately, both the nucleons may 
enter the target, form a compound assembly where they excite 
one, few or many nucleons and eventual] y a proton is released 
again. The life time of the compound state may be expected to 
be small if fewer nucleons get excited by the incoming deuteron 
while that of a state which is complicated may be expected to 
be relatively long. One would therefore expect resonances with 
varying widths on the basis of uncertainty principle. In the 
event that the probability of the simplest process, where the 
proton does not enter the nucleus at all, is large compared to 
others, a treatment purely on the basis of direct reactions 
might give satisfactory results. In general, however, both pro- 
cesses will contriuute to the observed proton yield. This is 
the basis for adding the amplitudes as calculated on the basis 
of two respective mechanisms. In general the two amplitudes 
have to be added coherently when the wave-packet representing 
the incident particle is much longer than the expected time 
delay caused by the direct and compound nucleus reaction 
mechanisms. A resonance of width aEq is observed with the 
beam of energy resolution AEb, when 


AEg is smaller than . 
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The corresponding uncertainty At„ in the arrival of the beam 
particle is larger than the time delay At^^ which is of the 
order of nuclear life time. Hence the dimension of the beam 
packet IS large compared to the distortion produced at its edges 
by the finite life time of the compound nucleus. It is expec- 
ted therefore that the compound -nucleus amplitude will add 
coherently and produce interference effects The experimental 
evidence of this situation is provided by elastic scattering 
experiments where one finds interference effects between poten- 
tial and resonant scattering in the excitation curve. 


Direct Reactions 

In the direct reaction theory the Hamiltonian of the 
system is considered in the form^^ 


H = t ...(2.9a) 

= H(^ + Tjj + Vq ... (2.9b) 

— Hg + Tg + Vg ... (2 . 9c) 


where Hq is the sura of the internal Hamiltonians of the two 
clusters a and A, T^ the relative kinetic energy of the two 
clusters and is the sum of the potential energies of all 
particles m the cluster a due to each particle in cluster A. 
The form (2.9) of the Hamiltonian is not suitable for conside- 
ration of bound wave functions but is suitable for the 
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consideration of unbound wave functions. Let$g be the eigen 
function of H 3 so that, 

Hg ...(2,10) 

The Schrodinger equation for the total wave function is 

H = E ... (2.11) 


taking the scalar product with we have 
($g, H = E (Og, 4'^+^) 


. .( 2 . 12 ) 


where the integration is only over the internal variables of 
Substituting Eq. (2.9c) for H into (2 12) we get 

(Tg+Sg-E) ($ 3 , = -($p, Vg 'J'^'*'^) ...(2.13) 

Since the scalar product ( ) involves integration over 

only the internal variables o-f iJg , it defines a function of 
r-g , which IS denoted by 

(«e' 'h > 


Therefore, we have 


(Tj +^- E) 5^ = -(»j, Vg fW) 


. . . (2.14) 


This equation has solution 




iKq . r„ 




2irh 


ft iKg(fg-r^) , , __ 

^ I “ # I P P P 




. . . (2.15) 
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where we have 


E -e 


Kg)' 

e" W7~ 


Taking the asymptotic term of eqn (2 15) , the amplitude 






27rh" 


( e 


iKg . rg 


, Vg ) 


. . (2.16) 


The formula represented by Eq (2.16) is exact. It can 
be recast into a different form which is useful for the dis- 
cussion of applications We add the term V- ( 4>- , to both 

P P 

sides of Eq. (2.13)/ where Vg is an arbitrary function of r^ . 
This leads to the relation 


(Tg + Vg + Eg - E) = -($g, (Vg-Vg)'F^'^^) . .(2.17) 

The solution of this eqn. is obtained by using the Green's 
function technique as above. This leads to the following 
expression for the amplitude 


where 


(Kc,/ Kg) 



= _ J!? ( $ Y (Y - V ) 

2 TT n 

is a solution of the homogeneous Eq. 
(Tg + Vg + Eg -E) Xg = O 


. . . (2.18) 


. . . (2.19) 


Both forms given by Eq. (2.17) and (2,19) are equivalent. For 
applying either of the tv;o forms a knowledge of the total 
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wavefunction 'i' is needed. In the direct reaction theory 
the bound state component of is neglected From amongst 

the unbound comnonents often only the component in the incident 
scautering channel is retained. Use of more than one unbound 
component leads to a more complicated coupled channel calcula- 
tion, . The calculation further requires the use of the channel 
potential Vg . Simplifying assumptions are usually made for 
calculations . 

In the Plane Wave Born Approximation (PWBA) the total 
wave function is replaced only by the incident plane wave 

e in the incident channel a The scattered part of 

ot 

the incident wave 3s therefore further omitted. In the case 
of (d,p) reactions the approximation to Vg is made as follows. 
We have 


V. = V + V „ 
p pn pN 


where V is the proton neutron interaction and V is the 
pn pN 

total interaction of the outgoing proton with all the nucleons 
of the target. In the pure stripping mechanism model of Butler 
it IS assumed that the proton does not enter the nucleus at all 
It IS therefore assumed that the contribution due to the inter- 
action could be neglected in evaluating the integral (2.16) 
In this approximation the amplitude becomes 
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Actg - - 




-ikg.rp 

2 / *8 ® '^np *“= 


2irh 


— / < 4 1 <)i ^ e 

2 J t'g I 


-ik . r ik ^ . r T 

P ^ V 4), e ^ dr dr 
np d II P 


... (2 . 20 ) 


where 4>p ^-nd cji^ are expressed in terras of the proton and the 
deuteron coordinaces respectively. The integration in the 
scalar product 1 over the internal coordinates of the 

target nucleus A and it has the form 






s (r .,) 
1 nN 


. . . (2 20a) 


where s is essentially the coefficient of fractional parentage 

and T), (r .J is a normalised function of the neutron coordinate 
1 nN 

typical of the angular momentum 1 carried by the neutron bound 
to the target N 

Further the interaction V in Eo. (2.2o) can be written 

np 

in terms of deuteron binding energy by writing the Schrodinger 

equation for the deuteron. Finally, changing the variables from 

r , r to the relative coordinates r and r , we get 
n p pn nN ^ 


O y -iq.r 

2 N pNrr i /- \ _ npy 


A= s + q dr^p e 




iQ.r 


nN 


. . . ( 2 . 21 ) 
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where q and Q are essentially the change of momenta of the 
proton and the nucleus v/ith some scale factors. In order 
to evaluate the integrals in (2.21) the exponential factors 
are expanded in terms of spherical Bessel functions as follows 

e^5*^ = A-kI (JJg) ...(2.22) 

X ]i 

Since involves spherical harmonics of order zero and rji 
that of order 1^ the integration over angles leads to the 
appearance of Bessel functions and respectively. For 
the evaluation of integral over r^^ often a cut off radius R of 
the order of nuclear radius is used. The angular dependence 
of the amplitude is then primarily given by 

A 'U Jq (QR) ... (2.23) 

The angular distribution of the proton is therefore 
determined by the value of 1 and the momentum transfer Q, 
which depends upon the angle of emission. The radius R is 
usually used as a parameter within limits. The Bessel func- 
tion has an oscillatory behaviour with rapid damping. The main 
peak of !j^(QR)i is compared with the experimental mam peak 
in (d,p) reactions. The peak position in many (d,p) reactions 
IS predicted to a good approximation by this theory but not 
in all cases. The agresnent away from the peak is however 
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not found satisfactory. For prediction of absolute cross- 
sections estimate of the coefficient s is obtained from 
structure consideratn on and using this the predicted cross- 
section IS found to be higher than the observed cross-section 
by a large factor 

In the Distorted Wave Bom Approximation (DWBA) the 
total wave function is replaced by the incident channel 

wave function including the scattered part. To obtain this 
incident channel wave- function , scattering data in the 
incident channel is used. The parameters of the optical model 
potential are first determined using scattering data and 
these are then used to calculate the incident channel wave- 
function. In the expression (2.18) for the amplitude^ the func 
tion Vg is any arbitrary function of the channel radius r^ and 
the function x„ is then an eigen function of the two body prob- 

p 

lem with the potential . In the DWBA the function is 
chosen equal to optical model potential in channel g which 
leads to the observed scattering data in the channel g . For 
the (d,p) reactions we then have 


(Vg-Vg) = V 


pn 


+ V 


pN 


- V 


opt 


. . . (2.24) 


It IS then assumed that the last two terms cancel each other 
to a good app^oxxnation, so that the remaining interaction is 


then V. 
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The interaction V is the sum of all interactions 

pN 

that the individual nucleons of the target exert on the out- 
going proton. It IS not clear that this sum can be represen- 
ted by a function of the channel radius alone. Such reduction 
of a many body interaction to an equivalent two body inter- 
action IS an approximation. However, the reduction of the 
operator (V^-Vg) in (2.18) to by a proper choice of Vg is 
a much better approximation than neglecting as is done in 
the PWBA. Thus the DWbA method involves a better approxima- 
tion to well as to the operator in the matrix eleraent- 

The amplitude in Eq. (2.2o) is then given by 




^0 


2Trfi 


2 •' 3 *3 pn a 


. (2.25) 


where the functions X^are the solutions of 


^2, 2 

r fi k 

“ 2Vq 


T - VJ X = O 


6^ ^ = 


. (2.26) 


with the asymptotic forms, 


(+) ik.r \ e 

► e + f (e, a) 


ikr 




-+ e^^*^ + f (‘n’-B, TT + cti) 


-ikr 


. (2.27) 


The optical potentials Vg are of the form, 


^3 = ^c + 
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where U is the Coulomb potential in the channel and U is 
^ IJ 

the optical potential. This potential is of the form 

Uj, = + iWj 

and various forms for the real and the imaginary part are 
commonly used. 

The amplitude (2.25) can be written in the form 

*06 = -J / aii / ai-g Xj 'S|Vp^|a> xy ...(2.28) 

where J is the Jacobian of transformation to the relative co- 
ordinates, and <3 I a> includes integration over the inter- 

nal co-ordirates of the target. 


<3 Iv 1 a> 

pn 


= < I V i 4). d> > 

B ' pn^ ^a 


= 1 V 


s n(r -^)v d) (r ) 
' nN pn a np 


This factor, which depends upon r^^ and r^^ can be expressed 
in terms of the channel radii r^ and r^ . The integral (2.28) 
IS then a SIX dimensional integral Often the interaction 
IS assumed to be of zero range and is replaced by a delta 
function 6 • This effectively leads to the 
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simplification/ 

- 

-nr = ^ 

So that the integral (2 2 8) reduces to three dimensions 

In evaluating the integral (2.28), the matrix element 
<3lVp^la> IS expanded into a series of multipoles each of 
wnich corresponds to transfer to rhe target nucleus/ a definite 
angular momentum j composed of an orbital part 1 and a spin 
part S. We have then 


j = 1 + S = Jg - 

and in the absence of the spin-flip, we have 

Often only one value of j/ 1 and S is important or allowed in 
a transition. The distorted waves x also expanded in terms 

of partial waves, 


4Tr 


' 17 ^LM ’‘l 


we then have 


1 - Lq - Lg 


The cross-section, which is proportional to the square of the 
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amplitude, can then be written in the form 


22 


da _ ^ 


Ib P a 

' Isi ' DW 


d<!- 2J^+1 isj 723^15 


- . . (2.29a) 


where 


a Isj 
DW 


4ir 




M. 


B 


M 


'A 


mmamb 


1.0x10 T r 

2s +1 { 21 +l)i 


. . . (2.29b) 


where B, is the measure of strength and is proportional to the 

IS j 


. , nimam-r^ 

spectroscopic factor S 2 .SJ* quantity ° rs given by 


’isj 


I 




Xj l 

a b a £) 




(J^j m^-mra i ^a’^a^ (2L^+1) (L^loo| L^o) v{2s^+l) (2i+l) 


\/(2J-|^+l) (2Lg^+l) 


r^b ^b ^b 
1 s j 

^a ®a '^a 


d 0 .^) 


M. 


A 


M 


i~Tr~~ / Xt (h, ~ r ) f-i (r ) Xt (hr)) 

' Mg •' a Lb r Mj^ a Is j a a a 


. . . (2.30) 


where m = mo-m. +m, - m and d (k, ) is the rotational 

BA D a ribma~m b 

matrix. The symbol ( | ) represents the Clebsch-Gordan coeffi- 


cients and f-, (r ) is the form factor depending upon the 
xsj a 
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choice of the neutron bound state wave function and Ib^ 

' Isj 

IS proportional to the interaction strength and the 
spectroscopic factor . 

In the case of (d,p) reaction the formula (2-29b) 
reduces to 


do 

dn 


1.53 


21 ^ 4-1 



'Isi 

(2j+ 1) 


Isj 


. . . (2.31) . 


2 3 

where Hulthen function has been assumed for the deuteron 
wave function. Normally only one value of Isj is required 
and therefore (2.31) simplifies to the form 


dp 

df2 


2J^+1 


s « 
^Isi 

(2j + 1) 


Compound Nucleus Reactions 

A discussion of formation of compound nucleus would 
really need the consideration of the time development of the 
non- stationary state produced in nuclear reactions. The state 
at t= - “develops into a bound state at the instant of com- 
pound nucleus formation which subsequently decays providing 
flux in unbound channel in an exponential manner. In the 
compound nucleus model it is assumed that the compo- 

nent of the stationary state wave function predominates within 
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the nuclear volume to provide the boundary conditions at the 
nuclear surface and the stationary unbound wave functions 
predominate just outside the nuclear volume. The coupling 
between the two is provided by the boundary conditions. The 
compound nucleus has a definite spin and parity so that 

we have 


loi + - Ig + Sb + Jg 

^a \ (-1)^ = '^3 “ '^c ...(2.31a) 

It IS usual to introduce channel spins such that 

= \ + + ..{2.31b) 

As a conseguence of restrictions (2.31) the values of angular 
momenta in the channels a and g are restricted by 

la = (Jc + ' * *'^c " ®cxl 

l0 = (Jc + • I'^c “ ^sl ...(2.32) 

Further, because of parity conservation the permitted values 
of 1 q and Ig are either all odd or all even. The angular 
dependence of the amplitude in channel 0 therefore, has 

the form 

^iTiax , 

A„d = E a y 

1 ( even or odd ) ^ ’ 

m 


(2.33) 
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The cross-section is proportional to the square of the ampli- 
tude and therefore has terms which are of even order. The 
dependence on the azimuthal angle vanishes after summation 
over the magnetic quantum numbers and the cross-section can 
be expressed in the form 


dP 

d n 


L 

max 

I 

L 

even 


P- (cos 


©) 


. . . (2.34) 


after expansion in terms of Legendre polynomials. The angular 
distribution^ therefore, are symmetric about 90° in the centre 
of mass coordinates. If more than one level of the compound 
nucleus of different panties contribute the summation (2.33) 
will be over even as well as odd values of 1 and the angular 
distribution will then no longer be symmetric about 90° CM. In 
the case of contribution from two or more levels of the com- 
pound nucleus the angular distribution changes rapidly with the 
bombarding energy, since the relative contribution from diffe- 
rent levels changes. 

The angular dependence of the amplitude A , 

s m s ra , 
s s' 

from initial channel s m_ to final channel s' mi is given by 

S o 

2 a 

(Blatt and Biedenham) 


*^3 m^ s' m^ ' 
s s 


1 ' -1 J, 


- I (21+1)^ (Iso m^) 


1’ 


X d's' i 

... (2.35) 
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where T' 


sms' m , IS related to amplitude q ■ by 

s s' 


the relation 


{a 


Va H 


s m s ' m , 
s s 


(2irh ) 


2 V 2 ' ^s m s ’ m , 

/ en c* • 


Eq.2.35 is the detailed form of (2.33). Ihe summation is over 

all quantxim numoers except s, m , s', m , The quantity U is 

s s 

the partial reduced amplitude which is independent of angle 
but depends upon the system energy. If a single level of the 
compound nucleus is excited, there is no summation over J^. If 
a statistically large number of levels contribute it is 
expected that the interference terms due to their randomness 
may cancel and the cross-section will be relatively insensi- 
tive to 0. Also, angular distribution will be approximately 
symmetric about 90*^ in centre of mass system. The total cross- 
section IS obtained by summing over the final channel spins 
and averaging over the initial channel spins in case of un- 
polarized beam. 

The calculation of the quantity U is to be carried out 
as a problem of nuclear structure in the region of high exci- 
tation energy. This leads to the prediction of the position, 
angular momentum, parity and other properties of the levels of 
tre cc^pounc. nucleus. The compound-nucleus model developed 
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so far does rot attempt to calculate these from first princi- 
ples However/ assuming that these levels exist, their effect 
on the nuclear reaction process can be predicted with some 
simplifying assumptions If at excitation energy such a 

level exists, the probability of finding the incident particle 
a within the nuclear volume is expected to be high, since at 
this energy a bound state is predominantly excited. The condi- 
tion of high probability within the nuclear volume is satisfied 
when the slope of the radial wave function rR^(r) of the un- 
bound state IS zero at the nuclear surface. The energy at 
which a level of the compound nucleus is formed, is therefore 
also the energy at which the logarithmic derivative f^ at the 
nuclear surface is zero. The general quantum mechanical theory 
of scattering in the presence of absorption shows that the 
cross-section for scattering and for compound nucleus forma- 
tion peah for this boundary condition It is therefore possi- 
ble to understand appearance of resonances in terms of levels 
of compound nucleus by introducing boundary conditions at the 
nuclear surface- The behaviour of cross-section in the 
neighbourhood of resonance is found to be of typical dispersive 
type as given by the Breit-Wigner formula. This feature is a 
consequence of the fact that the logarithmic derivative varies 
reasonably smoothly with energy in the neighbourhood of its 
value zero. In that case it can be expanded in power series, 
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about the resonant energy and the retention of the linear term 

alone directly leads to Breit-Wigner formula which is of the 

^ 2 4 

form 

TJ 

2 ...(2.36) 

(E-E^)+ lA 

2 

The cross-section is proportional to U and for the 
shape of the resonance the parameter A is found to be equal 
to r/2 . where r is the full width of the resonance at half 
maximum. According to uncertainty principle the width is pro- 
portional to the rate of decay of the compound nucleus. The 
rate of decay can be further divided into partial rates of 
decay into various specific channels. We can thus define 
partial widths r , r etc corresponding to these partial rates 
such that r = £ Fg The energy dependence of the quantity U of 
Eq. (2.35) is 



1 ' s ' 


''^^nls ^nl's' 


(E-E„) + 2 


n 


. . . (2.37) 


Substituting Eq (2.37) in (2 . 35) , summing over all quantum 
numbers, the amplitude q has the form 



. . (2.38) 
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where n represents different resonance energies, represents 
dependence on parameters other than energy and the phase factor 
IS included in g . 

Con tri bution of Dire ct and Comp ound tiucleus Reaction s 

As mentioned earlier the actual problem of finding the 
asymptotic behaviour of the total wave function, which includes 
the bound and unbound parts with the coupling between the two, 
is treated in terms of models The amplitude is assumed to be 
made of a resonant part and non-resonant part. This has basis 
in terms of semi-classical consideration and is applicable 
when reasonably resolved Breat-Wigner type of resonances are 
found riding over general background in the measured cross- 
section. This IS particularly applicable in the case of reac- 
tions like the (d,p) reactions on light targets where well 
recognized stripping patterns of angular distributions are 
observed at a bombarding energy of about 4-5 Mev while at bom- 
barding energy of about 2-MeV the stripping patterns appear in 
a somewhat distorted form and resonances appear in the excita- 
tion curves . 

In such cases it is possible to assume that the ampli- 
tude q IS given by 

qi, (Q) 

a(E,©) = I — — - + (©) -..(2.39) 

{(E-EP+ 
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rif O 

where the non-resonant amplitude q may be considered to be 
energy insensitive in magnitude and phase over the width of 
the resonance, though it may depend upon 0 The phase of the 
resonant part of the amplitude varies with energy as deter- 
mined by Breit-VJigner formula. Due to this energy dependence 
it IS possible to empirically separate the non-resonant part 
of the cross-section from the resonant part by analysing the 
experimental excitation curve data at each anale It is then 
possible to reconstruct the angular distributions of the non- 
resonant mechanism alone and compare these with the available 
predictions of the direct reaction theory. 



CHAPTER III 

EXPERIMENT 


Accelerator 

The present work was done with a two-Mev Van-de- Craaff 
accelerator available at the Indian Institute of Technology, 
Kanpur. It is a suandard two-Mev machine of High Voltage 
Engineering Corporation, Burlington, Mass., U.S.A 

The accelerator was purchased without the gas drying 
system which was designed and assembled locally. The syston 
provides a fresh charge of nitrogen and carbon dioxide mixture 
with a dew point of -60°F in about ten hours which is then 
transported to the pressure tank at 350 psi to provide insu- 
lation of the high voltage dome. A water chilling plant was 
similarly locally assembled. 

The ion source is a standard radio frequency source 
which IS at present used with hydrogen and deuterium as source 
gases. The ions are focussed and accelerated down the beam 
tube to provide total currents upto about l50 A amperes with an 
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energ]^' resolution of about lo KeV. The analysing and switch- 
ing magnet for beam handling was under construction while this 
worlv was in progress. Therefore, a laboratory magnet provid- 
ing about lo hilo- gauss of magnetic field was used to separate 
the molecular beam from the atomic beam. 

The voltage on the dome was measured by means of a 

generating voltmeter (GVM) fixed inside the pressure tank. It 

consists of a rotor which alternately shields and reveals a 

stator to the electric field produced by the charge of the 

dome. An AC signal is therefore generated on the stator which 

is amplified, rectified and measured. The output DC signal is 

proportional to dome voltage. The absolute calibration of the 

'1 

beam energy was done by locating the Li(p,n) threshold expec- 
ted at 1,88 MeV and the ^"^Al(p,Y) resonance expected at 0.992 
MeV. The neutrons were detected with a BF^ counter embedded 
in paraffin and Y“3rays were detected with a Nal crystal of 
dimension 5"x4" mounted on photomultiplier. The final dome 
voltage was read on a digital voltmeter. The average fluctu- 
ation in the dome voltage during the period of data collection 
was +3 KeV . 

Target Chamber 

A scattering chamber with a diameter of 10" was designed 
for charged particle experiments and is shown in Figure 1. A 
system of two slits was used to collimate the beam to a narrow 
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spot of two millimeters in diameter and a third slit with a 
larger aperture was placed nearer the target as an anti- scatte- 
ring baffle. The target was mounted on a brass holder which 
could be raised, lowered or rotated without breaking the vacuum 
by adjusting the clamps on the top of the scattering chamber 
The centering of the target could be ensured by visual checks 
of the fluorescence through the perspex window provided for the 
purpose. The scattering chamber had the provision of holding 
a number of solid state detectors at any desirable radial dis- 
tance from the target and their angular position could be read 
and varied without breaking the vacuum with a knob at the 
bottom of the chamoer. Arrangament for a monitoring detector 
at a fixed angle was also made. An electrically isulated Para- 
day cup was provided to collect the transmitted beam. This is 
fitted with a secondary electron suppressing ring which was 
kept at -150 volts to prevent leakage of charge through secon- 
dary electron emission from the Paraday cup. The output of the 
Paraday cup was fed to an Eloor current integrator which 
measures the current as well as the integrated charge. The 
chamber was provided with an evacuating port and a vacuum of 
lo” mm was maintained using the evacuating system of the 


machine. 
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Detector 

For the detection of char gai particles surface barrier 

solid state detectors were used. They had a depletion layer 

sufficiently thick to stop 6 MeV protons and 24 MeV alpha 

particles. These were obtained from the Oak Ridge Technical 

Corporation (ORIFC) , Oak Ridge, Term , u.s.A and some were 

obtained from the Bhabha Atomic Research Centre, Trombay, 

Bombay. Some detectors made in this laboratory were also used 

for preliminary studies. The detector signals were fed to ORTEC 

model lol and 102 low noise preamplifiers and model 201 and 202 

ORTEC amplifiers The output of the amplifier was analysed 

using a 512-chennel multi-channel analyzer (MCA) purchased from 

Nuclear Data Corporation. The resolution of the detection sys- 

241 

tern was checked with the Am alpha source and was found to 
be 0.8% for the 5.48 MeV alpha particles. It was improved to 
0.7% v;hen only the central portion of the detectors was used 
with a detector collimation system. A slight effect of degra- 
dation of resolution, due to neutron irradiation, was observed 
during the course of (d,p) experiment. The MCA itself was 
calibrated using a linear pulser provided in the ORTEC ampli- 
fier system 

Target 

12 9 

The targets of C and Be were made by vacuum evapo- 
ration technique. The system was assembled in this laboratory 
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and IS shown in Figure 2. The electrode system couplings and 
the liquid nitrogen trap were fahricated locally. It was possi- 
ble to obtain a vacuum of 10 mm in the evaporation dome and 
a current of about hundred amperes could be used for evapora- 
tion of the material. 

For the preparation of carbon targets, two carbon elec- 
trodes (99.9% purity) were fitted to the live posts. The tip 
of one was tapered and that of the other was a plane surface. 

The tips were brought in touch with a rotary seal arranganent 
to strike an arc. A tantalum boat was fitted to two other 
live posts and a small quantity of sodium chloride was placed 
in the boat. Glass slides etched with hydrofluoric acid and 
cleaned with acetone were held at about 15 cm above the boat 
for deposition. A clean glass dome was placed above the whole 
system and the dome was evacuated. The sodium chloride was 
first evaporated by passing current through the boat and the 
glass slides were thus coated with thin-layer of sodium chloride. 
An^arc was then struck between the carbon electrodes which led 
to the deposition of a carbon film over the sodium chloride 
film. The slides were then removed and slowly immersed in dis- 
tilled water which dissolved sodium chloride. The carbon films 
floating on the surface of water \rere then picked up on rectan- 
gular brass mounts with circular holes in them. 
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For preparing the beryllium targets carbon films were 
first prepared as oacking and the glass slides were replaced 
by these mounted films for evaporation of beryllium. Beryllium 
(99.99% purity) was evaporated using tungstun filament instead 
of a tantalum boat, since beryllium forms an amalgum with tan- 
talum. The beryllium targets with carbon backing were used in 
the experiment and the additional peaks in the reaction spec- 
trum arising out of beryllium reactions could be easily detec- 
ted since the carbon peaks were already studied. The beryllium 
peaks of interest were well separated from the carbon peaks due 
to the difference in the Q values. Attempts to prepare a self 
supporting beryllium film did not succeed since the beryllium 
flaked-off v/hen attempt was made to peel off the film from the 
glass slides. 

The relative measurement of the thickness of the carbon 
targets used was obtained from their respective yields of 
(d,p^) reaction. The absolute thickness was found by weighing 
a relatively thick film in a microbalance and also by studying 
the energy loss of 1.90 MeV protons by studying the shift in 

n 

the threshold of Li (p/n) reaction as shown in Figure 3. The 
two methods gave consistent results to within 10%. The carbon 
targets used had a thickness between 3-16 KeV for 1.9 MeV protons. 
In the case of beiryllium target, the thickness of the backing 
carbon film was first determined and that of beryllium film 
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was obtained by measuring the energy loss through the two 
films together. Beryllium targets used had a thicbness of 
3 to 7 KeV for 1.9 MeV protons. To improve an overall accu- 
racy^ thick-target thickness was determined. 

Measu rem ent s 

For alignment of the beam, the Faraday cup was replaced by 
a pyrex disc and the position of the scattering chamber was 
adjusted to ensure that the beam spot was at the centre of the 
pyrex disc. The targeb was then lowered and the passage of the 
beam through its centre was ensured by checking visually the 
bright spot on the target appearing due to beam bombardment. 

The pyrex disc was then replaced by the Faraday cup. The zero 
degree position for the detector was checked by observing the 
yield of the protons scattered from carbon 12 at the same posi- 
tive and negative angles for several angular positions. 

A typical spectrum at 1600 KeV deuteron energy and 
®lab ” 113° beryllium on carbon backing is shown in Figure 4. 
All resolved peaks in the spectrum could be identified from the 
known Q values for the respective reactions For reaction angles 
smaller than 70° lab , the elastic scattering cross-section in- 
creased rapidly as compared to reaction cross-sections. This 
increased the dead time of the equipment to an undesirable 
extent. Therefore, a nickel film was placed in front of the 
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detector to just stop the elastically scattered deuterons but 
not the energetic protons from (d,p) reactions- With nickel 
foil in front of the detector, data upto 15°Lab. could be taken 
in the forward direction. A typical spectrum with nickel foil 
in front of the detector is shown in Figure 5- To identify 
peaks in this case, loss of energy in Ni foil for each reaction 
product had to be calculated separately. 

12 

The excitation function for c(d,p ) reaction was 

obtained at nine laboratory angles in steps of 10 Ke/ in the 

deuteron energy range of 800 to 1850 KeV, shown in Figures 

l2 

6-9. The excitation curve for c(d,d) channel, in the same 
range, is shown in Figure 10. Further, for consistency check, 
angular distributions in the angular range 15*^ to 160° labo- 
ratory angles in steps of 10° were taken at seventeen values 
of bombarding energy in the range 800 to 1800 KeV. These are 
shown in Figures 11 to 13. In the low energy range of 500-800 KeV 
no resonances were observed. Excitation curve in this range at 
the three laboratory angles and angular distributions at two 
bombarding energies are shown in Figures lo and 14. 

On ^Be target data on (d/P^) , (d,t), and(d,a^) 

could be obtained. Since these reactions have a higher Q value 
than the cCd/p^) reaction, the corresponding peaks could be 
studied inspite of carbon backing at most of the angles. In 
the energy range of 800 to 1800 KeV excitation functions for 
%e{d,p^) reactions were taken at nine angles. These are 
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shown in Figures 6 , 15 and 16. For consistency checTc angular 
distribution in the angular range 15° to 160° laboratory 
angles in steps of 10° were taken at eighteen values of bom- 
barding energy in the range 500 to 2000 KeV as shown in 

g 

Figures 17 to 20. In addition, excitation fxinctions for Be 

(d,p^) reaction in the range 800 to 1800 KeV and excitation 
9 

function for Be (d,cc^) in the same bombarding energy range 
are shown in Figures 21 and 22. Excitation curve for Be {d, t^) 
in the range 800-1800 KeV is sho-'/n in Figure 23. 

It IS estimated that the cross-sections have been mea- 
sured with an energy resolution of 0 , 5 %. The relative error 
in the cross-sections is about 2% and the error in absolute 
cross-section is about 12% in carbon and 15% for beryllium. 

In the present analysis the accuracy in relative cross-section 
IS of primary importance. 
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CHAPTER IV 


analysis of data 

12 

The excitation functions of C {d,p^) reaction leading 

13 

to the ground state of c at nine angles is shown in Figures 

6 to 9. The excitation function of G (cl,d) reaction at two 

backward angles is shown in Figure 10. The striking feature 

in this data is the appearance of correlated peaks at all 

angles in (d/P^) reaction as well as in the (d yd) channel. For 

example there is a rather well separated peak appearing at 

about 1450 KeV in all these excitation curves. Interference 

effects are obviously present and these lead to the difference 

of details of the excitation curves at different angles. Thus 

the 1450 KeV resonance is sharpest at 140° in the (dyd) channel 

and is v;eak at 80° in the (d,p) channel. Other peaks v/hieh are 

less marked stand out to different extent at different angles. 

Thus the peaks in the vicinity of 930 KeV, 1310 KeV, 1770 KeV : 

and a broad hump around 1200 KeV can be noticedy which is con- 

17 

sistent with the previous data. Tne data indicates the 
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presence of a resonant process with a level density which is 
not loo hinh ^his ig expected for the relatively light com- 
pound nucleus n at an excitation energv of about 11 5 MeV 
The data obtainel m the past has been analysed by assuming 
the exiKjtence of compound nucleus levels at appropriate exci- 
tation energies 

The angular distribution of the ^^c(d p^) reaction in 
the energy range covered is shown in Figures 11 to 13 The 
noticeable feature of most of these angular distriburions is a 
maximum at about 35° and a strong backward peaking The maxi- 
mvm at about 35° is typical of the direct deuteron stripoing 
reaction while the backwa d peaking is not The angular dis- 
tributions on and off 1450 KeV resonance are shown in Figure ?4 

The angular distributions do not change much in their general 

! t- 

shape however the backward cross-section is large at the peav 
energy The dashed curves shown are those predicted on the 
basis of PJane Wavv. Butler theory with 1^^= 1 for the captured 
neutron The fits obtained for the 35° maximum are satisfac- 
tory However the theory grossly underestimates the relative 
cross-section at backward angles and at the ninimum of the 
angular distribution The available data at higher energies 
shows that the backward peaking in the angular distribution 

16 18 

decreases until more familiar stripping patterns are observed 
(shown in Figure 25) 
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The excitation function of ®Be(d,p^) reaction leading 
to the ground state of ^^Be at nine~angles is shown in Figures 
6/ 15 and 16 The excitation function of protons from the 

Be(d/P^) reaction leading to the first excited stete of ^^Be 
IS shown in Figure 21 The excitation functions of ^Be(d 
and "^Be(d t) leading to the ground state of "^Li and ^Be respec- 
tively are shown m Figures 22 and 23 Correlated weah reso- 
nances or antiresonances are observed in the excitation curves 
mentioned above at deuteron energies of 950 KeV, 1080 KeV 
1200 KeVy 1300 KeVy 1450 FeV 1620 KeV and 1760 KeV The reso- 
nances in all these cases are weak and shov/ that compound nucl^s 
effect IS snail as compared to those in the case of C (dyp^) 
reaction The excitation energy of the compound nucleus 

nil T A 

Be + d s= ■‘B IS about 17 MeV which is higher than in the F 
compound state Many channels are open in this case for the 
decay of the compound nucleus and therefore contribution to 
any one channel as expected to be small On tne other hand 
stripping process/ if present, will contribute predominantly 
to (dyp) and (dyn) channels Earlier Canavan has analysed 
^Be(dyp ) data in this energy range assuming that compound 
nucleus process alone is responsible for the reaction mecha- 
nism Later groups concluded that the reaction was taking 

place mainly through direct mechanism and tried to fit data 

20 

using the stripping theory 
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The angular «3istritoutions of the Be(d,p^) reaction 
-it about 2 ileV is shovm in Figures 17 to 20 In the forward 
angles a maxirnum at aoout 35*^ is observed at most of the ener- 
gies Also an equally strong backward peaking throughout is 
observed The peal at 35° is tynical of strinoing reaction 
mechanisn The angular distribution changes a little as the 
energy is changed However/ the relative forward to bad ward 
cross-section increases with the energy A typical angular 
distribution at 1310 KeV, 1445 HeV, and 1620 KeV is shown in 
Figure 26 Dashed curves are Plane Wave Butler curves with the 
neutron captured in the 'p* orbit The fits are good in the 
forward angle but predictions of relative cross-section remains 
poor in the backward angles The angular distribution data at 
higher energies * (Figure 27) shows considerable decrease 
in the relative backward cross-section and a pure stripping 
type of pattern with a small backward cross-section is observed 
at these energies 

In view of these observations it was felt that the 
possibility of coherent contribution from both the resonant 
and non— resonant processes should be investigated# although / as 
was done in the past/ it is possible to fit the angular distri- 
butions with Legendre polynomials and explain them in terms of 
contribution from a number of compound nucleus levels exclu- 
ding aprior the ross-ble ro’^— resorant contribution 
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In this connection the angular distributions obtained 

were fitted by us in tems of Legendre polynomials. It was 

found that a satisfactory value of in the least square 

fitting could be obtained by using both even and odd terms 

urto a maximum order of L, given by L = 6 in the case of 

max 

12 o 

C (d,p^) and 7 in the case of "BeCd/p^) We have 


where 1 is the smaller of the two values 1^ and Ig. In the 

12 q 

case of C we therefote have Iq 3 and in the case of Be 

la 4. To exolain the observed distributions on the basis 
of the compound nucleus theory rather hign partial waves in 
the incident beam have to be invoiced. On the other hand it 
IS exoected that sach high partial waves would lead to stri- 
pping mode rether than compound nucleus formation. 


Separation of iv- Resonant Amplitude 

The separation of compound nucleus part from excitation 

16 27 

curve has been earlier suggested by Bonner et al. and Yoshida 
The separation of non-resonant part from the excitation curves 
IS done on the basis of its smooth energy dependence. The 
differential cross-section in terms of the amplitudes 

3 s ' 

IS given by the relation 


da I 

dn I sra^s’m , 
s ^ 



I q 1 

' ^sm s m , 


2 


... (4.1) 
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The amplitudes q„ ^ ^ ^ is split into a non-reso~ 

o fii S ill t 

s s 

nant (NR) and a resonant (R) part At a given angle 9 / we 
have 

qCE^) = + q^(E„) ...(4.2) 

where the channel quantum numbers have been dropped for the 
convenience of v/riting the expressions The resonant ampli- 
tude has the following energy dependence 

D 1 6 

C ^ e ^ 

q^(Ea) = I . . (4.3) 

n (E -E ) + 1 r„/2 

CL n 


where the summation is over all the contributing resonances 
R 

and ' IS independent of energy The quantity e is the 

tlli 

relative phase due to n resonance. 


When only one level of tne compound nucleus is consi- 
dered and in Eg. (4.3) contribution from single resonance is 
considered, the cross-section is given by 


dg 

dg 


sm s 
s 






(E -E )+ ir^/2 
o o o 


2 


...(4.4) 


The amplitude q^^^" is assumed to be constant over the width of 
the resonance. Equation (4.4) therefore becomes 
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df'" 
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HR 


s m s ' p , 
s s 


r, 2 I 


a 


o 




^ j ^FR R j 2 cos 5^(E^-E^)+ 2T^/^ sin 6 ^ 

(F^-s^)’-+ rj/ 4 

^2 ...(4.5) 

\/here factor — has been absorbed in the ampljtudes. To 

evaluate the cross-section ^ , the Eq. (4.5) is sumried over 

the final channels and averaged over the initial channels. 

The relative contribution of each term in Ea (4.5) does not 

change by this averaging. Carrying out prooer summation and 

re-defining the constants an (4,5) after terms of similar 

deoendence are grouped together, we get. 


^ ^ D 

ao E„ 




+ r^V4 Ei, ((Ei-E^)hr^V4) 


(4.6) 


The constants D, , and are independent of energy 
Ejj. The first term D/s^gives the value of non-resonant part 
of the cross-section. The Eq. (4.6) can be re-written as 


Y = 


I 


P f 
2 2 


(E, 




^ ^3 ^3 


(E, 






. . . (4.7) 


The Eq. (4,7) becomes a basis for separation of non-resonant 
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contribution from total cross-section, when contribution from 
cnly one resonance is considered The constants P^, and 

are treated as parameters in least squares fitting of this 
expression (4 7) 


If cor tr Lbutioii fron two levels of the compound nucleus 
IS considered, Sq (4.3) is re— wrir.t'^n as. 


g(E ) + 


C ^ e ° 


P 

=1 


(E^,-E^) + ir/2 


. . (4.8) 


The cross-section --p: , 

dn 's m s ' m , 
s s 


, is given by 
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.. (4,9) 
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Eq. (4.9) IS further simplified by collecting terms of similar 
energy dependence. After summing over the final channels and 
averaging over the initial channels, the differential cross- 
section IS given by 


^ ^ D_ 

^0. E„((E„-E^)h 


h tJ^/A E„ ( V4) 


+ 




+ 


E, 


{ (E„-E^) +r^ /4) ( (E^ -E^) ^+r^^/4) 


C(E„-E^)^ + T^^/4) C{E„-e/+ r^V 4 ) 


((e^-Eq)^ + r^V 4 ) C(e^-e^)2 + r^V 4 ) 


. (4.10) 


The constaits E, <^ 0 ' , C^. have been re-defined 

to include all other constants. The relative phase 6^, and 
are also treated as constants for a given 6 and meraed in these 
constants. The function (4.10) can be rewritten as 



.. (4,11) 
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Expressions (4.7) and (4.11) could noM be. used for least 
SQuares fitting to the experimental excitation curves taken 
at various angles. 

T/i/hile in principle it was possible to fit the excita- 
tion curves in the entire energy reaion studied by considering 
the contribution of each resonance everywhere, the number of 
terms in the fitting function and the number of parameters 
become large. Instead , a few regions of the excitation curve 
around relatively isolated resonances were studied. In the 
first step the contributions from the resonance under consi- 
deration and the contribution from the non-resonant process 
were considered The formula (4.7) was therefore used for 
least souares fitting in this case, ■'•''ext the contribution from 
one neighbouring resonance was added and the least square 
fitting was done this tine with formula (4.11). This was 
extended to consider the contribution to a single resonance 

by the non -r esonax' c process and by the two neighbouring reso- 

do 

nances The aim at each step was to obtain the value of 
given by the first term of the fitting functions at each 
angle. in the cases studied it was found that the non-reso- 
nant cross-section as obtained by using function (4.7) changed 

when function (4.11) was used instead, but changed little 

dO' 

thereafter. Therefore for theextraction of ^ it was 
found sufficient to consider the contribution to the mam 
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resonance from the non-resonant process and one neighbouring 
resonance. The resonant energies and the widths of resonances 
were varied along with the coefficients in ■’he fitting func- 
tion, to obtain the best fit. By this method it ^/as found 
that the excitation function in the neighbourhood of 1450 FeV 
in the case of C(d,p^) could be accounted for by resonance 
of width lo KeV located at 1445 KeV, a non-resonant contribu- 
tion and the contribution from a resonance of width 40 keV 
located at 1310 KeV. The fits obtained in this region of exci- 
taticn curve are shown in Figure 28. The non-resonant contribu- 
tion at 1445 KeV was determined at each angle at which excita- 
tion curves were obtained. It was then possible to plot the 
angular distribution of the non-resonant part alone at 1445 FeV 
bombarding energy. This extracted non-resonant cross-section 
was then compared with the DWBA theory of stripping. This 
procedure was repeated at the 930 KeV resonance whose width 
was determined as 110 FeV A contribution from neighbour reso- 
nance located at 1160 KeV was considered Width of this reso- 
nance v/as detenraned to be 190 KeV The non-resonant cross- 
section was also extracted at 1770 KeV with 140 KeV width- 
The extracted angular distributions due to non-resonant process 
are shown in Figure 29. The non-resonant contribution accounts 
for about 20% of the observed cross-section. 
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In Be the non- resonant cross-section contribution was 
determined at 920 KeV, 1310 KeV, 1450 KeV, 1620 KeV and 1770 
KeV, where widths of resonances were identified as lOO KeV, 

50 KeV/ 80 KeV, 50 KeV and 120 FeV. The angular distributions 
due to non-resonant contributions are shown in Figure 30. The 
non-resonant cross-section in this case accounts for about 
90^0 of the observed cross-section 

DWBA Calculations 

The extracted angular distribution due to the non-reso- 
nant process was compared witn the direct reaction theory of 
the stripping process The DWBA calculations were carried 
out with program DWUCK ana the computer BESb-16 at Bhabha 
Atomic Research Centre, Trombay, Bombay 

The distorted wave program DWCK calculates the diffe- 
rential cV^Oss -sect ion as a function of angle 0 for a given 

ItS 1 

energy of the incident particle. The quantity as 

given by Eg. (2.29b) is calculated for a given value of 
orbital angular momentum transfer 1, spin transfer s and 
total angular momentum transfer j . The distroted wave 
for partial wave L is calculated by solving the following 
differential equation niomerically * 

(Vg )] X^(r> = O 

dr'^ r -fi 


. . . ( 4 . 12 ) 
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where 


Vc 




SO 


and the potenrial l^^(r)is qiven by 

6fiX) 




'p " '^p' -dx: 


. . . (4.13) 


where Vl^ is real part of the potential and is the imagi- 


nary part. The function f(X) is given by 

>- 

f(X) = (1 + 


X = (- 


r-R A 
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1/3 


-) 


a = Diffuseness parameter 
A = Atomic Number of the nucleus. 


The potential gives Coulomb interaction due to charged 
sphere of radius and charge and is given by 

Z Z e^ 

U = ■ - — r > R 

or c 


■J "a 
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(3- r <P^ 

R ^ 
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. . . (4,14) 


Tins xntsicsctioii U IS jjx'opox'txon.^-l to spin oxtoxt xn.tsx'sctxori 

SO 


L,S and IS given by 


vsoR 1 r - 

^SO = h 4572 ? — P— 


dr 


... (4,15) 
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where VSOR comes from the Thomas Spin Orbit term and has a 
value of 25 for protons and neutrons. The wave function 
satisfies th^ boundary condition in the asymptotic 
region where j 0 and also inside the nucleus where 

X^(0) ^ 0 The are ccnputed by numerically integra- 

ting the differential ecfUation (^.12) and thus values of 
Xj^(r) for a set of discrete values of r is obtained in the 
entrance and the exit channel with proper choice of potential 
parameters in the channel concerned. In the case of stripping 
reaction where a particle is captured in the outer orbit, the 
form factor is obtained by calculating the bound state wave- 
f unction for the captured particle such that it reproduces 
correct binding energy. TIius for stripping 

where normalized radial wave-function of the 

captured particle. After the calculation of form factor 

f^ (r), and distorted wave Xt ^ integral in Eq. (2.30) is 
Isj ' ‘ ' L 

obtained by doing numerical integration. Finally program 
calculates the quantity as given by Eq (2.29b) by 

inserting proper values of Clebsch Gordan Coefficients and 
spins of the system. 

There is more than one choice for the form of potential 
U(r) in Eq. (4.12) The commonly used form is given in 
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Eq. (4.13)« Here real potential is Wood-‘3axon type and imagi- 
nary potential is derivative of this function. 

The calculation of integral in Eg. (2.30) is done using 
zero-range approximarion. The DWUCK program also has facility 
for taking into account correction due to finite range effect 
and IS calculated in the following v/ay. 

The overlap of wave function and 4’^ given by 

r(x) = / d 

2 

and the fourier transform G(K ), is defined as 

iK.X 

G(K ) ~ fe D(X) dX 

The zero range approximation is given by 
Dq = G(0) 

and finite range correction parameter FNRHG is evaluated by 

2 2 

taking first term in the expansion of G(K ) about K — O. that 
IS 

_2 1 9 g(K^) 

R — ' 5 T 

G(K^) 3(K^) j.2^ ^ 

Thus in the reaction where 

a = b + c and B = A + c 
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the finite range correction is computed by the quantity 


W(r) = { 1 + 


2 \ W 




M. 


. (4.16) 


Potential Uj^(r)in differential Eg. (4 12) is a local 
potential. The program DWUCK takes into account correction 
due to non-locality. The effect of the non-locality is taken 
into account by multiplying a wave function for an equivalent 
local potential Uj^(r)by the factor^® 


f (r) 


(1 


7Z^ 

4 ^2 


U, 


N 


(r))^ 


. .( 4 . 17 ) 


un the channel under consideration. Here 3 is range of 
non-locality and is 0.85 for protons, 0.54 for deuteron and 
0.2 to 0.3 for ^He. 

The effect of finite-range correction and non-locality 
in general results in reduction of the cross-section, but the 
over all change in the shape of angular distribution is not 
expected. 

In calculating the DWBA cross -sections wath DWUCK, the 
principal input consists of the optical model parameters in 
both the incoming and oi-rgoing channels. "Ilie values of optical 
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parameters which provides satisfactory fit at high energy of 
2 5 2 9 3 

about 5 MeV ' were used as intial values in fitting the 
extracted curve. These values are given in Table 1. 


Table 1 


Target 

Incident 

Particle 

Real 

Potential 

Im 8 gina 2 ry Potential 
( Derivative type) 

^^I ^I ^I 



^R 

^Be 

D 

-160 

0.9 

0.9 

-12 

2.1 

0.5 

r^= 1.3 

P 

-49 

1.25 

0.65 

-7 

1.25 

0.47 

1 = C 

n 

-69 

1.60 

0.61 

-4.0 

1.41 

0. 65 

1.3 

P 

-51 

1.25 

0.65 

-7.0 

1.25 

0.47 


In the first place, since the deuteron bombarding 
energy is lower than the Coulomb barrier, it was thought that 
the nuclear interaction of the optical potential may not be 
necessary and Coulomb distortions alone would be sufficient 
for the calculations. In this case it may be noted that 
according to Eg (2.18) the fully distorted waves, including 
nuclear di stortions/must be used in the final channel 6 
although the wave function could be tried as the incident 
wave with the Coulomb distortion alone. This is necessary if 
Vg-Vp IS to be equated to Preliminary DWBA calculations 
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with only the Coulomb interaction in the incident channel 
were carried out. The calculated angular distributions were 
somewhat similar to the Butler curves v/ith higher cross-sec- 
tions predicted at backward angles. One such curve for 
c 

■^Be(d,p^) at 1445 KeV is shown in Figure 31 The fit for 
the backward angles is poor. A more serious discrepancy is 
however m the absolute cross-section where the predicted 
value IS about 15 times the observed value It was therefore 
concluded that a nuclear distortion of the incident wave can 
not be neglected. 

12 

In the DWBA calculations for C (d,p^) the optical 
parameters were varied about their initial values . Further 
calculations were made with and without spin orbit term wirh 
finite range correction and possible correction due to non- 
local effects were also investigated. The inclusion of spin- 
orbit term in rhe optical potential increased the predicted 
back angle cross-section. The inclusion of corrections due 
to finite range and non— local effect did not alter the shape 
of the angular distribution significantly but reduced the 
overall cross-section (Figure 32). For ^Be(d,p^) reactions 
only finite range correction was included. Figures 33 to 36 
show DV^BA fits to ^^C(d/P^) with only spin orbit term. Figures 

37 to 41 show DWBA fits to ^BeCd.p^) The optical parameters 

12 9 

which gave the best fit in the case of C and Be are given 
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in Table 2. 


Table 2 


Target 

Incident 

Particle 

Real Potential 

Imaginary Potential 
(Derivative type) 





^I 

^I 

^Be 


D 

-160 

0.9 

0.9 

-12 

2 1 

0,47 

r = 

1.3 








c 


P 

-49.5 

1.25 

0.65 

-7.0 

1.25 

0.65 



D 

-69 

1. 6 

0.61 

— 4 » 0 

1 — I 

0.65 


1.3 

P 

. 

in 

1 

1.45 

0.65 

-5.4 

1.22 

0-65 


It may be noted that the optical parameters given in 
Table 2 are close to those given in Table 1 for higher energy 
data. 

Expression (2.31) for the direct cross-section for 

Is 1 

(d,p^) reaction involves the spectroscopic factor s . 

A comparison of the DV©A calculation v/ith the experimental 
results provides a value of spectroscopic factor s which 
IS to be compared with that obtained frcan nuclear structure 
calculations.^^ The values obtained from the rWBA analysis 
of the data and those obtained from nuclear structure calcula- 
tions are shown in Table 3. 
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Table 3 


Energy (KeV) 


Experimental 


Calcd. 

920 

1310 1445 1620 

1770 

^Be(d,p^) 

2.2 

1.9 2.0 1.9 

1.7 

2.36 

^^C(d,p^) 

3.0 

1.2 1.1 

1.2 

0 613 


It may be noted that the DWBA analysis of the unextrac- 
ted experimental curve would have led to a spectroscopic factor 

which would be larger than the theoretically expected values 

12 

by a factor of about 8 for the case of C(d,p^). 

We may summarise the results of the analysis as follows 

12 

The excitation functions in C(d,p^) reaction show resonances 
which appear at various angles and also in the (d,d) channel. 
The angular distributions indicate the possible existence of 
a non-resonant process also. It is possible to extract the 
non-resonant cross-section from the observed excitation curve 
on the basis of Imown energy variations of the resonant ampli- 
tude (Breit-Wigner foimnula) and the slow energy variation of 
the non-resonant amplitude. The non-resonant contribution is 
about 20%. The extracted angular distributions have been 
compared with the DWBA. calculations made with optical para- 
meters close to those used at higher energies. The agreement 
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in the shape and absolute value is satisfactory* The results 
on Be(d,p^) are similar The resonances are^however, weaher 
and the direct contribution is found to be about 909' • This 
IS understandable since due to higher excitation of the con- 
pound nucleus , more chdnnels are open and a corresponding 
smaller percentage of the compound nucleus contribution is in 
the proton channel Although the fit in the backward direction 

n 

in the case of "^Be is not satisfactory, the forward cross- 
section in the angular dependence and magnitude is well 
reproduced . 



CHAPTER V 


CONCLUSION 


The aim of this work was to understand the reaction 

mechanism of (d,p) reactions on light targets at sub-Coulomb 

bombarding energies in the range of 0.8-2 0 MeV. The targets 
g 12 

used were Be and C The features of the stripping mecha- 
nism such as the small angle peak for 1=1 capture, which are 
easily identified at about 4 MeV bombarding energies, were 
still present thouah not as distinct as in the high energy 
data, indicating the possible presence of the stripping mode 
at low energies The excitation functions however, showed 
resonances typical of compound-nucleus-formation mode. Assuming 
a coherent resonant and non-resonant contribution it was possi- 
ble to separate the non-resonant cross-section and compare it 
wibh that predicted on the basis of DWBA calculations of the 
stripping mode. Prom the analysis the following conclusions 


could be drawn 
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(1) The aor roach of accounting for the asymptotic 
behaviour of the unbound wave-functions / v/hich are coupled 
to the bound wave functions/ in terms of coherent contribu- 
tions obtained from direct reaction model and compound -nucleus 
model (based on physical considerations) leads to satisfactory 
agreement v/ith the data in the reactions studied. 

(2) The contribution of direct reactions to the cross- 

section in (d/p ) reaction is about 20‘/a and rn ^Be (d/p ) 

o o 

it IS about 9C% in the energy range studied. 

(3) Even at sub-Coulomb bombarding energies used, it 
was found that the distortion of the incident wave is impor- 
tant and the Coulomb distortion alone does not lead to satis- 
factory results. 

(4) At higher bombarding energies of about 4 MeV the 
stripping mode predominates and the optical model parameters 
for the DWBA calculations have been well determined. Essen- 
tially the same parameters can be used to fit the angular dis- 
tributions at the energies studied provided that the resonant 
part is extracted out. 

The effect of finite- range correction and non-local 
correction in the ^^C reaction is to reduce the cross- 

section in magnitude by about 1C9^ without much change in the 
shape of the angular distribution. 



Ill 


(6) Resonances were located in the compound nucleus 

1 4 17a 
N at excitation energies consistent v;ith earlier results. 

The widths of the resonances which we obtained are however 

somewhat on the lower side in most cases. This difference 

15 to be expected since we have included non-i'esonant contri- 
bution in the analysis. The spin and parity assignment made 
previously by neglecting the non-resonant contribution needs 
re- examination . 
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